Loss and revival of phase coherence in a Bose-Einstein condensate 
moving through an optical lattice 
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We investigate the phase coherence of a trapped Bose-Einstein condensate that undergoes a 
dynamical superfluid-insulator transition in the presence of a one-dimensional optical lattice. We 
study the evolution of the condensate after a sudden displacement of the harmonic trapping potential 
by solving the Gross-Pitaevskii equation, and comparing the results with the prediction of two 
effective ID models. We show that, owing to the 3D nature of the system, the breakdown of the 
superfluid current above a critical displacement is not associated to a sharp transition, but there 
exists a range of displacements for which the condensate can recover a certain degree of coherence. 
We also discuss the implications on the interference pattern after the ballistic expansion as measured 
in recent experiments at LENS. 

PACS numbers: 03.75.Kk, 05.75.Lm 



I. INTRODUCTION 

The possibility of manipulating Bose-Einstein conden- 
sates (BECs) in periodic potentials has provided the op- 
portunity to investigate a wide range of phenomena, ex- 
ploring a very stimulatingfield which combines atomic 
and solid state physics 0, B H II S IS S 13 ■ 

In particular, the dynamical behavior and coherence 
properties of BECs loaded in optical lattices have been 
the subject of an extensive experimental 0, IE IE 13 
and theoretical work jE HE HE ■ Recently it has been 
demonstrated that a trapped BEC moving through a one- 
dimensional optical lattice created by a laser standing 
wave can realize a Josephson junction array sustaining 
an oscillating atomic current 5] . The condensate is set in 
motion by a sudden displacement of the trapping poten- 
tial, and then, when the center-of-mass velocity reaches 
a critical value, the system undergoes a transition from 
the superfluid regime to an insulator regime, character- 
ized by a localization of the condensate in the trapping 
potential jj. 

This phenomenon, which is accompanied by a loss of 
coherence, is triggered by the onset of a dynamical insta- 
bility, as theoretically demonstrated in 9,'To', TT|. These 
studies rely on the analysis of the Bogoliubov spectrum 
of simplified one-dimensional models, and predict a sharp 
transition from the superfluid to the insulator regime, 
with a complete loss of the system coherence. Such a be- 
havior has been confirmed in the deep insulator regime, 
that is for large initial displacements, by the direct solu- 
tion of the three dimensional Gross-Pitaevskii equation 
(GPE-3D) 13. 

In this work we show that, owing to the 3D character 



of the system, a non-trivial phenomenology takes place in 
the intermediate regime where the condensate is able to 
retain a certain degree of coherence. The role played by 
the dimensionality of the system is investigated by com- 
paring explicitly the solution of the GPE-3D with those 
of two one-dimensional models. We also discuss how the 
phase coherence of the system influences the interference 
pattern after the ballistic expansion, showing that the in- 
terference peaks which characterize the superfluid phase 
01 are not completely destroyed in the insulator regime, 
though they may have a reduced visibility. The onset 
of decoherence processes is signaled also by additional 
structures ("fringes") which appear in the central peak, 
as recently observed at LENS j7tlll|. 

The paper is organized as follows. In the next section 
we introduce the GPE-3D and the two one-dimensional 
models considered (GPE-ID and NPSE). Then, in sec- 
tion lllll we present the results by discussing the evolution 
of the phase coherence, the Fourier power spectrum of the 
system, the center-of-mass dynamics, and eventually the 
effect of the free expansion of the system. The conclu- 
sions are drawn in section Hvl 



II. GENERAL FORMULATION 

Let us consider a Bose-Einstein condensate trapped in 
a cylindrically symmetric harmonic potential superim- 
posed to an optical lattice. The dynamics of the system 
is described by the three-dimensional Gross-Pitaevskii 
equation (GPE-3D) ^ 
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where N is the number of condensed atoms, g — 
A-nh^a/m the coupling strength, m the atomic mass and 
a the inter-atomic scattering length. The external poten- 
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with 



FIG. 1: Schematic of the system setup. 



tial V{x) — Vho{^) + Vp{z) is the sum of the harmonic 
trapping potential 



(2) 



and of the periodic potential generated by the optical 
lattice along the axial direction z 



Vp{z) — S Er COS^ 



27rz 



(3) 



A being the wavelength of the laser, Er = /2mX^ the 
recoil energy of an atom absorbing one lattice photon, 
d = A/2 the distance between two adjacent minima (lat- 
tice sites), and s a dimensionless parameter controlling 
the intensity of the lattice. 

To model the LENS experiment 01 , we consider a *^Rb 
condensate with = 5-10* atoms, and the following 
parameters characterizing the external potential: = 
2-K ■ 92 Hz, = 27r • 9 Hz, A = 795 nm, and s = 5. 

The ground state of the system in the combined har- 
monic plus periodic potential is found by mapping the 
wave function on a discretized grid 0| and using a stan- 
dard imaginary time evolution 'l5l. Then, at t = 0, the 
harmonic potential is shifted by Az and the condensate 
is set to evolve in the combined potentials as sketched in 
Fig. n To solve the time-dependent GPE-3D we use a 
split-step method which combines a FFT evolution in the 
axial direction p?7l fl8i and a Crank-Nicholson algorithm 
for the radial one [l9l| . 

In order to investigate the role played by the transverse 
dimension we also compare the full three-dimensional be- 
havior of the GPE-3D with the solutions of two one- 
dimensional effective models that account for the axial 
dynamics along the lattice direction. 

The first model considered is the ID Gross-Pitaevskii 
equation (GPE-ID) given by 
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where giD is obtained by a suitable renormalization of the 
3D coupling constant g. In particular, by requiring the 
invariance of the Thomas- Fermi (TF) chemical potential 
(that is the invariance of axial size of the condensate in 
the TF hmit) |23|, we get 
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where = ^JhJ(jnuj^ and aj^ — y'?i/(m are the 
characteristic oscillator lengths. 

The other model considered here is described by the 
Non-Polynomial Schrodinger equation (NPSE) [21| ob- 
tained from the GPE-3D by means of a factorization of 
the condensate wavefunction in the product of a gaus- 
sian radial component of width a{z,i), and of an axial 
wavefunction t/;(z,t) that satisfies the differential equa- 
tion (NPSE) 
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coupled with an algebraic equation for the radial width 

<j{z, t) = ax ^l + 2a,7V|i/;(z,i)|2 . (9) 

Owing to this partial coupling between axial and radial 
degrees of freedom, the NPSE is expected, with respect 
to the GPE-ID, to provide a more accurate description 
of the actual ground-state and dynamics of the system 
at least in the coherent regime, as discussed in 21, 22]. 

III. RESULTS AND DISCUSSION 

In this section we investigate the phase coherence and 
the center-of-mass dynamics of the condensate during the 
trapped evolution through the optical lattice. To discuss 
the degree of coherence of the system we also show the 
evolution of the momentum distribution of the conden- 
sate and the expected signatures after the free expansion 
of the system. 

A. Coherence 

In order to characterize the degree of coherence of the 
condensate inside the optical lattice we define the quan- 
tity x(^) the squared modulus of the correlation be- 
tween values of the condensate wave function evaluated 
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FIG. 2: Coherence of the condensate between nearest- 
neighbor sites of the optical lattice, obtained from the so- 
lution of the GPE-3D, for several initial displacements. After 
the onset of the dynamical instability (for Az ~ 40 /im) the 
system undergoes a sudden loss of coherence, that is partially 
recovered for later times. 



at points corresponding to the distance between nearest- 
neighbor sites of the lattice ^3 



(fx **(r, z; t) ^'(r, z + d;t) 



(10) 



In the foUowing we wiU refer to x '"^ "coherence" : x — ^ 
means that the system is fully coherent, whereas x = 
indicates a complete decoherence and the loss of the 
superfluid properties. 

In Fig. 121 we show the evolution of the coherence for 
several initial displacement, up to 80 ms. The picture 
shows that above a critical displacement, after an initial 
coherent evolution, the system undergoes a sudden loss 
of coherence, that is partially recovered in the subsequent 
evolution. We have verified that for longer times (up to 
200 ms) the value of x(i) does not grow anymore and 
remains close to the "saturation value" at 80 ms. 

We note that even though the system may become 
completely incoherent (% — 0), coherence is preserved 
in each sub-condensate (in each lattice site), and this is 
sufficient to justify our description of the system in terms 
of solutions of a differential equation. 

As shown in Fig. lathis capability of the system to re- 
gain a certain degree of coherence is due to its 3D nature. 
In fact, both the NPSE and the GPE-ID predict a com- 
plete dephasing after the onset of the instability, charac- 
terized by an almost vanishing coherence. The failure of 
these ID effective models in the insulator regime is due 
to the fact that at the onset of the instability most of the 
energy initially associated to the coherent evolution of 
the system is absorbed by the modes responsible for the 
dynamical instability, which grow exponentially in time 
destroying the phase coherence of the condensate 0, 0| . 
However, while in the ID case all the energy is transferred 
to these modes, in the 3D case the interplay between ra- 
dial and axial degrees of freedom, coupled through the 
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FIG. 3: Coherence of the condensate between nearest- 
neighbor sites of the optical lattice for a displacement = 
40 ^m. When the system enter the insulator regime the ID 
models (GPE-ID, empty circles; NPSE, triangles) predict a 
complete loss of coherence, while the actual (3D) behavior of 
the system is characterized by a partial recover of coherence 
(GPE-3D, filled circles). 



nonlinear term, allows for a partial transfer of energy to 
coherent modes of the condensate, acting, as a matter of 
fact, as a restraint to the decoherence process [l^. Note 
in this respect that the partial coupling between axial 
and radial degrees of freedom in the NPSE is not suffi- 
cient to account for the rephasing mechanism, and this 
is likely due to the simple gaussian factorization of the 
wave-function which produces a coupling just between 
radial and axial densities [2lj |. losing any information on 
the phase dynamics psf. 



B. Center-of-mass dynamics 

Let us now discuss the center-of-mass evolution along 
the lattice direction, by considering the axial coordinate 



Zcm{t) = / (fx z\^{r,z-t)\'' 



(11) 



and its velocity 
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In Fig. H we show the evolution of Vcm{t) during the 
trapped dynamics, for several initial displacements (see 
also Fig. 121). This picture points out clearly the existence 
of a critical velocity independent of the initial displace- 
ment (here |ucrit| — 1.25 mm/s, in nice agreement with 
the experimental value reported in Q), beyond which the 
system cannot sustain a coherent oscillation. This is due 
to the presence of axial modes with imaginary frequency 
that grow exponentially in time and are responsible for 
the dynamical instability of the syst em in a certain range 
of condensate quasimomenta |9|, Il0| . 



4 



-0.1 



-0.2 







' 




- Az=28 \\.m 
Az=36 \\m 
Az=40 
Az=44fim 

- ■ Az=52 \\.m 

- Az=60 


x;^.^ Vc'--'^"-' /■ 





20 



40 
t (ms) 



60 



80 



FIG. 4: Center-of-mass velocity of the condensate (in units of 
the Bragg velocity vb ~ hn/{md)) obtained from the solution 
of the GPE-3D, for several initial displacements. This picture 
clearly shows that the decoherence process is characterized by 
the achievement of a critical velocity (see also Fig. 1^. 
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FIG. 5: Center-of-mass dynamics for three values of the ini- 
tial displacement {Az — 36 fim continuous line, Az — 40 fim 
dashed-dotted, Az — 60 /im dashed). Contrarily to what 
happens in the simplified ID models (GPE-ID, empty cir- 
cles; NPSE, triangles), in the insulator regime the condensate 
center-of-mass continues a slow evolution toward the center 
of the trap (GPE-3D, filled circles). 



In Fig. Ulwe compare the corresponding center-of-mass 
evolution with the prediction of the one-dimensional 
models, for three initial displacements. We notice that 
the GPE-3D predicts a slow motion toward the center 
of the trapping potential, as observed in the experiment 
0,^3' whereas the ID models predict a sudden localiza- 
tion of the system as a consequence of the complete loss 
of coherence. 

Note also that for small initial displacements, that is in 
the regime of coherent oscillations, the predictions of the 
NPSE are in nice agreement with those obtained by the 
full GPE-3D, while the GPE-ID deviates from the ex- 
pected behavior after about 40 ms of evolution (we have 
observed in a similar fashion that the GPE-ID slightly 



underestimates the critical displacement). This effect 
may be due to the chosen renormalization of the cou- 
pling constant gi (see Eq. ©) that can overestimate 
non-linear effects [22| responsible for the existence of dy- 
namically unstable modes. However, even though one 
tunes gi in order to better reproduce the actual critical 
displacement, this would not change the inadequateness 
of the model to describe the behavior of the system in 
the insulator regime, neither the fact that in any case the 
NPSE represents a more reliable effective model (with re- 
spect of the "simple" GPE-ID) to describe the dynamics 
of the system in the superfluid regime pH . 



C. Momentum evolution 

Let us now consider the evolution of the axial momen- 
tum density distribution (power spectrum), which pro- 
vides a deeper understanding of the system behavior and 
of the dephasing mechanisms. In Figs. 1^1 and [T] we show 
the axial power spectrum of the condensate, defined as 



P{p,) = 2n rdr\^{r,p,)\^ 



(13) 



(the tilde indicates the Fourier transform along z), in the 
case of an initial displacement Az = 40/xm, at subsequent 
evolution times. In particular, in Fig. I^we consider three 
configurations when the condensate is still in the super- 
fluid regime, respectively for t = 0ms (a), 15ms (b) and 
30 ms (c) . In this case the power spectrum is character- 
ized by sharp peaks localized at Pz{t) = ±2npB + q(t) 
(n = 0, ±1, . . ., here only the zero and first order are vis- 
ible) ™, ps = mvB = fiTr/d being the Bragg momentum 
and q{t) the condensate quasimomentum. In fact, when 
the lattice intensity is sufficiently high (i.e. we are in the 
tight binding regime), the condensate in a Bloch state of 
quasimomentum q can be written as |23j (neglecting for 
simplicity the presence of the harmonic potential which 
breaks the periodicity of the lattice p^) 



\E'(r, z) 



k 



"^'"^ip{r, z + kd) 



(14) 



where ip(r, z + kd) are wave functions localized at each 
lattice site (labeled by the index k). It is straightforward 
to show that in the coherent regime the momentum dis- 
tribution is characterized by sharp peaks whose weight 
is modulated by the axial Fourier transform of the wave 
function at each lattice site 



sin^{Ni{p,-q)d/2n)^ 
sir?{{p,~q)d/2fi) 



2t: rdr\ip{p,,r)\^ (15) 



Ni being the number of occupied lattice sites 

Contrarily, the behavior of the system in the insula- 
tor regime is completely different, as shown in Fig. [3 for 
t = 45ms (a), 60ms (b), and 80ms (c). Indeed, after 
the onset of the dynamical instability, the momentum 
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FIG. 6; Axial power spectrum of the condensate in the su- 
perfluid regime for t = 0ms (a), 15ms (b) and 30ms (c). The 
presence of sharp peaks spaced apart by multiples of ±2ps is 
a signature of the system coherence Q. 



FIG. 7: Axial power spectrum of the condensate in the insu- 
lator regime for t = 45ms (a), 60ms (b), and 80ms (c). 



distribution initially spreads out losing any signature of 
coherence (Fig. [TJa)). Afterwards, the system tries to re- 
arrange itself in a coherent fashion as signaled by the ap- 
pearance of structures localized in correspondence of the 
initial peaks. Even though these "peaks" have a rather 
large spread (see Fig. Efc)), the relative population with 
respect to the central one is of the same order of that in 
the full coherent regime (Fig. El^a)). Notice also that the 
momentum distribution is centered not exactly in = 
since the condensate is slightly moving toward the trap 
center. 

For comparison, in Fig. |S| we show the axial power 
spectrum as obtained from the NPSE, after an evolu- 
tion of 80 ms. In this case there is no evidence of any 
localized structure (contrarily to the case of Fig. [7Jc)) 
as further confirmation of the fact that these effective 
one-dimensional models cannot account for an accurate 




p^ (units of q^) 



FIG. 8: Axial power spectrum of the NPSE solution after an 
evolution of 80 ms. The distribution is completely spread out 
indicating a complete loss of coherence. 
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description of the system behavior in the insulator regime 
(similar considerations hold for the GPE-ID). 



D. Free expansion 

In this section we discuss the expected behavior of the 
system after a ballistic expansion. When the system is in 
the superfluid regime the internal coherence of the con- 
densate shows up in a clear interference pattern charac- 
terized by lateral peaks that move outwards at velocities 
V — ±2pB/m 4|. Contrarily, in the insulator regime we 
expect the overall shape of the interference pattern and 
also the expansion of the central peak to be affected by 
the partial loss of coherence, as we will show in the fol- 
lowing. 

To investigate these aspects, instead of solving the full 
GPE-3D, we use a simplified model in order to avoid 
unnecessary heavy numerical computations p^ . In par- 
ticular, since we are mainly interested in the expanded 
axial profile, we neglect the contribution of the mean field 
interaction, and we assume a free expansion governed by 
the Schrodinger equation 



(a) 



2m 



«'(r,z;t) 



(16) 



The reason to use this approximation is threefold. 

(i) First of all, in our case most of the initial energy of 
the condensate (at the moment of the release from the 
trap) is associated to the fast density modulation, due 
to the presence of the optical lattice. Therefore, contrar- 
ily to "usual" case of a pure harmonic confinement, now 
the kinetic energy term predominates over the mean field 
one, E,nf ^ Ekin- 

(ii) Another reason is that we expect the mean-field 
interaction to affect mainly the radial expansion, which 
is however integrated out in our treatment. 

(iii) Finally, since our aim is to investigate the overall 
shape of the interference pattern and to point out the 
possible presence of lateral peaks (and/or other signa- 
tures of the degree of coherence of the system) rather 
than the exact size of the central peak, we can for this 
purpose neglect the contribution of the mean field term 

Then, by using Eq. Hlt)|) it is easy to show that the 
axial density distribution after an expansion time texp is 
given by (?i = 1) 



p{z;texp) = (fr\^{r,z;t. 



(17) 



^'0(7', Pz) being the axial Fourier transform of the initial 
wave function (at the time of the release from the trap) . 

As an example in Fig. |51we show the typical shape of 
the density distribution before (a) and after (b) an ex- 
pansion of about 30 ms, corresponding to a case similar 



I 





-400 -320 -240 



-160 -80 80 

z (|im) 



160 240 320 400 



FIG. 9: (a) Axial profile of the condensate in the trap, in the 
insulator regime (for a displacement Az = 40 /im and after 
70ms of trapped evolution), (b) Interference pattern (lateral 
peaks) and fringes (in the central peak) after a subsequent 
expansion of about 30 ms. The marked line in (b) is obtained 
by taking into account the finite resolution of the imaging 
apparatus as described in the text. 



to those in Figs. [ZJb)-(c). We see that both distribu- 
tions are characterized by a rather fragmented profile, 
although after the expansion (Fig. IHIb)) it is possible to 
identify some structures localized in correspondence of 
the first order peaks of the analogous expanded profile in 
the case of a fully coherent condensate (at z ~ 320 /im) 

i- 

Furthermore, to take into account the effect of the fi- 
nite resolution of the experimental imaging apparatus, in 
Fig. inib) we also show the convolution (marked line) of 
the axial density with a gaussian distribution of width a 
(here 2(t = 6 /^m, according to the typical experimental 
resolution). This procedure clearly evidences the persis- 
tence of lateral peaks and the appearance of fringes in the 
central one, whose presence is a signature of the partial 
loss of coherence (as the reduced visibility of the lateral 
peaks). 

Note that although the actual shape of the interference 
pattern depends on the initial conditions (as for example 
the trapped evolution time), the overall behavior shown 
in Fig. Elb) can be reproduced for a wide range of trap 
displacements and evolution/expansion times. We have 
also verified that the main features discussed in this work 
hold for lattice intensities in the range 2 < s < 10. 

Both the persistence of the lateral peaks and the ap- 
pearance of fringes in the central peak in the density dis- 
tribution of the expanded condensate have been recently 
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observed in the experiments at LENS 0, ^3 • 



IV. CONCLUSIONS 

We have studied the phase coherence of a Bose- 
Einstein condensate that undergoes a dynamical 
superfluid-insulator transition during the trapped evohi- 
tion in the presence of a one-dimensional optical lattice, 
as recently observed at LENS 0, [3| and discussed in 

From the comparison of the solution of the three- 
dimensional Gross-Pitaevskii equation with that of two 
effective ID models, we have demonstrated that the in- 
clusion of the transverse degrees of freedom is crucial to 
account for the actual behavior of the system, as observed 
in the experiments 0, . In particular we have shown 
that the breakdown of the superfluid current is not as- 
sociated to a sharp transition as predicted for the pure 
one-dimensional case, but there exists a range of param- 
eters for which the condensate can partially recover some 



coherence during the subsequent evolution. 

We have also shown that the degree of coherence of the 
system affects significantly the interference pattern after 
the ballistic expansion of the condensate, characterized 
by the persistence of lateral peaks (as a signature of a 
partial coherence) and by the appearance of fringes in 
the central peak (due to the dephasing of the system). 

These results open interesting questions about the pre- 
cise role played by the radial degrees of freedom on the 
excitation spectrum and on the decoherence mechanism. 
The investigation of these aspects requires the analysis 
of the 3D Bogoliubov spectrum of the system, and will 
be addressed in a future publication. 
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In the case of coherent excitations, and also for a free 
expansion, it is possible to use a more reliable ansatz as 
discussed in Ref. 0| 

This assumption is justified by the fact that the con- 
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